Entanglement versus relaxation and decoherence 
in a quantum algorithm for quantum chaos 
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We study analytically and numerically the behavior of the concurrence (a measure of the entan- 
glement of formation) of a pair of qubits in a quantum computer operating an efficient algorithm 
for quantum chaos. Our results show that in an ideal algorithm the entanglement decays exponen- 
tially with the diffusive relaxation rate induced by classical chaos. This decay reaches a residual 
level which drops exponentially with the number of qubits n q . Decoherence destroys the residual 
entanglement with a rate exponential in n q . 
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Enormous interest into quantum information and com- 
putation [jj] has generated serious efforts in characteriz- 
ing and understanding quantum entanglement which is 
considered as the ultimate origin of quantum power (see 
a recent review Q). A quantitative measure of the en- 
tanglement of formation, namely the concurrence C, was 
introduced and shown to be able to characterize an arbi- 
trary state of two qubits 0, |j] • Being closely related to 
the von Neumann entropy S of the reduced density ma- 
trix p of two aubits[3lj|. this quantity was recently found 
to have interesting applications to quantum phase transi- 
tions in interacting spin systems |5| . In parallel, the prop- 
erties of entanglement were investigated in a quantum 
model of coupled tops where it was shown that there ex- 
ists a typical value of entanglement which is determined 
by the chaotic behavior of the dynamics of the model 
@ and that the growth rate of entanglement of initially 
decoupled tops is increased by the underlying classical 
chaos On the same line, it was recently shown that 
contrary to intuition even a heat bath may create entan- 
glement between two qubits 0. 

All these studies |H, IS HI clearly demonstrated how 
rich entanglement properties can be in interacting quan- 
tum systems. However, in the context of quantum com- 
putation it is much more crucial to analyze the evolution 
of entanglement in a specific algorithm performing an op- 
erational task. Indeed, it is expected that the entangle- 
ment is very sensitive to noise and decoherence |!J Efl LLD] 
and the understanding of its behavior in an operating al- 
gorithm can lead to better strategies in the control of 
decoherence and imperfection effects. As for our knowl- 
edge such direct investigations have not been performed 
until now. Therefore, in this paper we study the behavior 
of the concurrence in an efficient algorithm for the quan- 
tum sawtooth map which has been proposed recently in 
|12|. The algorithm for this model has a number of im- 
portant advantages: all n q qubits are used in an optimal 
way and no ancillae are required, one map iteration in 



the Hilbert space of size N = 2 nq is performed in 0(n 2 ) 
quantum gates, and the algorithm is based on the quan- 
tum Fourier transform (QFT) which is one of the main 
elements of various quantum algorithms 0. This allows 
to simulate a complex dynamics in the regime of quantum 
chaos with a small number of qubits. Since the entangle- 
ment can be efficiently measured experimentally (see e.g. 
[l3|) the experimental observation of the concurrence be- 
havior discussed here can be realized on NMR 0, or 
ion-trap ^(| based quantum computers with about 6 - 
10 qubits and a few hundreds of gates. 

Contrary to the situation discussed in 0,0, our results 
show that in the exact quantum algorithm the underly- 
ing classical chaos leads to an exponential decrease of 
the concurrence C down to some residual level C which 
characterizes the global system coherence. On the other 
hand, the presence of noise in the quantum gates leads 
to a destruction of this coherence with a rate T growing 
exponentially with the number of qubits n q . This shows 
that entanglement can be very sensitive to decoherence. 

The dynamics of the classical sawtooth map ^3 is 
given by: 



= n — k 



dV(0) 
dO 



Tn mod 2n 



(1) 
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where V(9) = —9 2 /2, —pi < 9 < tt and the bars denote 
the variables after one iteration. After rescaling y = Tn 
and x — 8, it is clear that the dynamics depends only 
on the parameter K = kT. Due to the discontinuity 
in the derivative of V(9), the Kolmogorov-Arnold-Moser 
(KAM) theorem cannot be applied to the map and 
its dynamics becomes chaotic and diffusive for arbitrar- 
ily small values of the chaos parameter K > ^?|]- For 
K <C 1 the diffusion is governed by a nontrivial cantori 
regime which has been worked out in • I n this case the 
rescaled diffusion rate D (K) = (Ay) 2 /t w 1.2tt 2 K 2 - 5 /3 
is much smaller than the quasi-linear diffusion rate cor- 
responding to the random phase approximation D q i = 
tt 2 K 2 /3 (the latter becomes valid only at K > 1). The 
diffusion rate in n is D = Dq(K) /T 2 . 

The quantum sawtooth map [Hll9l.l20l.l2lL is the 
quantized version of the classical map, to which it corre- 
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FIG. 1: Dependence of the concurrence C on the dimen- 
sionless time t (the number of map iterations) for the map 
@ at K = 0.5, L = 4 and n q = 8, 12, 16 (curves from top 
to bottom respectively). The smooth curves show the fit 
C(t) = Aexp(— ■yt) + C of the relaxation to the asymptotic 
value C obtained on a larger time interval (t < 10 4 ). The 
inset shows C(t) on larger time scale. The initial state is 
(| 00) + |H))|0)/v2 where \<j>) is the uniform superposition of 
all but the two most significant qubits. Here and below the 
logarithms are decimal and all axis units are dimensionless. 

sponds in the limit k — > oo, T — > 0, and K = kT = const. 
One step in the quantum map is given by the unitary op- 
erator U acting on the wave function tp(9): 

^(0) = Uil>(6) = e jfce2/2 e- lT " 2/ V((9), (2) 

where n = —id/d6 (we set h = 1); when n q qubits 
are used for the discretization, so that the total num- 
ber of levels is N — 2 nq , we assume periodic bound- 
ary conditions on a torus in the phase representation 
(ip(9 + 27r) = ip(9)) and in the momentum representa- 
tion (ip(n + N) = ip(n)), as discussed in 0,|23|]. As a 
result, the phase 9 takes only N = 2 nq discrete equidis- 
tant values in the interval — 7r < 9 < n, and so does the 
momentum (n = 0, 1, N — 1, i.e. n + N is identified 
with n) . The sawtooth map approximately describes the 
dynamics in the stadium billiard [19j , the phenomenon of 
dynamical localization which is similar to the Anderson 
localization in disordered potentials and cantori induced 
localization [23, El • Thus this simple map describes 
a rich and complex dynamics and represents an interest- 
ing testing ground for efficient quantum computation. It 
is especially important to understand how this complex 
quantum dynamics will be affected by imperfections in 
realistic quantum computers. 

The numerical simulation of the map J5J is based on 
the quantum algorithm described in [121 ] and is imple- 
mented on the basis of a quantum computer language 
developed in |2J] which is well adapted for the experi- 
mental operation of quantum gates via classical computer 
software. In this way the dynamics of up to 20 qubits can 
be easily simulated on a laptop. 

To investigate the behavior of the concurrence in the 



quantum map J5J) we compute C for the two most signif- 
icant qubits which determine the first two binary digits 
ai : 2 in the expansion of momentum n: the reduced den- 
sity matrix p for this qubit pair is obtained by tracing out 
all other n q — 2 less significant qubits (the digits aj with 
3 < i < n q in the expansion of n — (0,10,20.3. ..a,i...a nq )). 
After that C is computed from p as described in 

Si- In 

this way we obtain the concurrence value C on a global 
scale of the whole system which is decomposed into 4 
equal parts with iV/4 quantum states in each of them. 
In addition we fix T = 2tt L/N in the regime of quantum 
resonance so that L gives the integer number of classical 
phase space cells embedded in the quantum torus of size 
N [T2I [23| (the classical dynamics is periodic in n with 
period 2ir/T). In the following we also take L to be a 
multiple of 4 to have an integer number of classical cells 
in the 4 parts of the partition in the momentum n. 

Typical examples of the dependence of C on the num- 
ber of map iterations t are shown in Fig.l. According 
to these data C(t) decays exponentially down to a resid- 
ual value C and, in the limit of large N, the decay rate 
7 becomes independent of N. It is natural to compare 
this rate with the rate of classical relaxation. Indeed, due 
to underlying classical chaos, the probability distribution 
fn = \^n\ 2 over n is described by the Fokker-Planck equa- 
tion df n /dt = Dd 2 f n /d 2 n/2 which gives the relaxation 
to equipartition with the rate 

7c - 2k 2 D/N 2 - D (K)/2L 2 . (3) 

The comparison between this classical value j c and the 
rate of concurrence decay 7 is given in Fig. 2. It clearly 
shows that the decay rate of C(t) is given by the classical 
rate: 7 = j c . It is important to stress that this relation 
remains valid also in the nontrivial cantori regime (K <C 
1) and that the quantum decay reproduces all oscillations 
of the classical diffusion (see inset in Fig. 2). 

The properties of the residual value of the concur- 
rence C are analyzed in Fig. 3. We will argue in the 
following that they can be understood in terms of the 
system conductance g ([2^, 26]); in view of this we ex- 
press C vs. g — 27 c / A = NDq(K)/L 2 where, up to 
a constant factor, A = \/N is the level spacing and 
27 c = D (K)/L 2 is the Thouless energy (see [2jJ for a 
recent review). In spite of strong fluctuations the data 
presented in Fig. 3 can be described by the global av- 
erage dependence C ~ ^ 1/VN. Indeed, for 
K = 0.5, L = 4 the system size varies by 3 orders of mag- 
nitude and the fit gives an algebraic decay with power 
a = 0.56 ± 0.02 being close to 1/2. We attribute the 
presence of strong fluctuations to the fact that the value 
C is averaged only over time but there is no averaging 
over parameters. Thus, from the point of view of disor- 
dered systems |27^ C represents only one value for one 
realization of disorder. 

We propose the following explanation of the results 
presented in Figs. 1-3. For a state \tp) like in Fig.l, we 
can write |i/>) = J2 ai a 2 l a i fl 2> \</>aia 2 ) where a lt2 = or 1. 
Then the value of the concurrence C is proportional to 
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FIG. 2: Dependence of the rescaled rate of the concurrence 
decay 7 = 2-yL 2 on the chaos parameter K for n q = 19, L — 16 
(triangles down); n 9 = 18, L = 8 (circles) and n q — 17, L — 4 
(triangles up). The solid curve gives the values of the diffu- 
sion rate Do(K) taken from Figs. 2, 3 of [17| showing that 7 
is determined by the classical relaxation rate The inset 
shows data on a larger scale with R — ~y/D q i (symbols) and 
R = D (K)/D q i (curve from 
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FIG. 3: Dependence of the residual value of the concurrence 
C on the system conductance g = NDo(K)/L 2 for a broad 
range of parameters: half filled circles show dependence on 
L = 4, 8, 12, 16, 20 for K = 0.5 and n q = 14, 15, 16; diamonds 
and triangles show the variation with K for n q — 14, L = 
16; n q = 15, L = 8 and n g = 16, L = 4. The filled circles 
connected by the dashed curve show the dependence on N 
for K = 0.5, L — 4. The solid line marks the slope X/-Jg. 

the difference of two scalar products, C ~ |Qi4 — Q23I, 

where Qu = 2 ^/| (0 O o |0n> I 2 and Q23 = 2y/\ (0oi l^io) I 2 - 
From this relation and the fact that the initial state is 
symmetrically distributed with respect to the transfor- 
mation n — > N — n, it follows that C is proportional to 
the difference \W\\ + Woo — W01 — W w \ where W ai a 2 is 
the total probability inside the part (0102). In the clas- 
sical limit this probability difference relaxes to zero with 
the classical relaxation rate 7 C , and that's why 7 = j c in 



agreement with the data of Fig. 2. 

The residual value C is determined by the quantum 
fluctuations of the previous difference of scalar products. 
In fact, due to the discretization of map J2J), the symme- 
try n — > N — n is broken and |0oo) becomes different from 
1 0n ) . Therefore in the scalar product Q14 (and Q23) the 
N/4 terms have random signs and thus Q14 oc 1/VN 
(each term is of the order of 1/N). In this estimate we 
assumed a summation over all N wave function compo- 
nents. However for finite values of the conductance g 
only the states inside the Thouless energy interval 27 c 
have a significant scalar overlap [27I l28j , and thus we can 
make a conjecture that N should be replaced by the ef- 
fective number of components which is of the order of 
N e g ~ 7c/ A ~ g. According to this, C ~ l/y 7 *? ln agree- 
ment with the data of Fig. 3. 

The existence of a residual level of concurrence for 
an ideal quantum algorithm reflects the fact that the 
global behavior of the whole system remains coherent. In 
fact, the Poincare theorem guarantees that for very large 
times the concurrence will have a revival close to the ini- 
tial value (however, this will happen on an exponentially 
large time scale). The situation becomes qualitatively 
different in presence of external decoherence represented 
by noisy gates. In our numerical simulations, noisy gates 
are modeled by unitary rotations by an angle randomly 
fluctuating in the interval (— e/2, e/2) around the perfect 
rotation angle. The presence of this external decoherence 
leads to a decrease of the residual value of C as illustrated 
in Fig. 4: the constant level is replaced by an exponential 
decay which gives C oc exp(— Ft). 

In order to obtain the dependence of T on the parame- 
ters we extracted it from the fit of the averaged ratio of C 
under a noisy evolution to its value in the ideal one. An 
example of such ratio and the corresponding fit is shown 
in Fig. 4. To suppress fluctuations we averaged over 20 
realizations of the noisy evolution. Moreover, the fit was 
restricted to the plateau regime, where the exact concur- 
rence is fluctuating around its residual value (the initial 
diffusive relaxation was excluded from the fit) . 

The results for T obtained in this way are presented 
in Fig. 5. Quite naturally we find that T oc e 2 , as it 
was also seen in other simulations of quantum algorithms 
with noisy gates (e.g. [29|). This scaling becomes better 
and better for large e values where T is larger. How- 
ever, more surprisingly there is an exponential growth of 
r with the number of qubits n q (T oc y/W). This result 
is very different from those obtained in 0, |2£j, where 
the time scale for fidelity and the decoherence rate for 
tunneling oscillations varied polynomially with n q . We 
see two possible reasons for the exponential sensitivity of 
the residual concurrence to decoherence. At first, in our 
case r is computed over a very large time interval, for 
which the quantum dynamics already reached its asymp- 
totic behavior (plateau for the residual concurrence); it 
is known that on very large times the eigenstates are ex- 
ponentially sensitive to imperfections due to the chaotic 
structure of the wave functions (see results and discus- 
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FIG. 4: Effects of decoherence on the residual concurrence: 
the two lower curves show the dependence C(t) for an ideal 
algorithm (top) and an algorithm with noisy gates at noise 
amplitude e — 0.003 (bottom). The time t is dimensionless 
(it is the number of map steps). In the latter case the average 
is done over 20 noise realizations. The curve in the upper 
part shows the ratio of C(t) at e = 0.003 to its value in an 
ideal algorithm, this ratio is averaged over a 100-kick moving 
window to reduce fluctuations. The dashed straight line shows 
a fit of the ratio to an exponential decay proportional to e~ Tt . 
Here n q = 12, K = 0.5, L = 4. 



sions in [3(| ). Another possible reason can be related to 
the fact that the residual value of the concurrence on the 
plateau is on its own exponentially small and maybe this 
is the reason why it becomes so sensitive to decoherence. 
Further investigation on the decoherence effects for the 
concurrence are required to understand in a better way 
this exponential sensitivity of C . 

In our studies we restricted ourselves to the investiga- 
tion of entanglement only between two qubits. The prob- 
lem of the characterizing the entanglement of a larger 
number of qubits represents an interesting challenge but 
at the same time it is much more complicated 0- How- 
ever even the relatively simple case of two qubits shows 
nontrivial links between concurrence and such interesting 
physical phenomena as quantum phase transitions [f| and 
statistical relaxation. 

In summary, our studies show that the decay of the 
concurrence in an operating quantum computer is deter- 
mined by the underlying relaxation rate of the classical 
dynamics. We show that the residual level of entangle- 
ment in an ideal algorithm scales as the inverse square 
root of the conductance of the system. This residual en- 
tanglement is destroyed by decoherence, whose effective 
rate grows exponentially with the number of qubits. 

This work was supported in part by the EC contracts 
RTN QTRANS and IST-FET EDIQIP and the NSA and 
ARDA under ARO contract No. DAAD19-01-1-0553. 
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FIG. 5: Dependence of the decoherence induced decay rate 
r of the residual concurrence on e 2 y/N for K = 0.5, L = 
4. Here the noise amplitude e changes from 0.001 to 0.01 
(10 equidistant values) for 7 < n q < 15. The data points 
(circles) are connected by lines for fixed values of e. The 
color intensity changes gradually from one chain to another 
to mark the variation of e (low/high intensity corresponds to 
small/large values of e). The straight line shows the averaged 
behavior T = 0.58e 2 V^- 



